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ABSTRACT 

,  In  order  to  determine  the  available  test  time  in  a  helium  driven 
shock  tunnel,  an  interface  detection  system  has  been  developed  to 
indicate  the  arrival  of  helium  driver  gas  at  the  test  section.  This  is 
done , by  exciting  the  gas  at  the  test  section  with  a  high  energy  elec¬ 
tron  beam  and  spectroscopically  analyzing  the  resulting  fluorescence  to 
determine  the  helium  concentration.  Relative  values  of  helium  and  air 

•  v:  ■ 

emission  intensities  are  calculated  and  compared  with  experimental  data 
from  both  static  and  transient  gas  mixtures. 
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CHAPTER  I 

TRANSIENT  GAS  CONCENTRATIONS  IN  SHOCK  TUNNELS 

'*1 


In  the  process  of  determining  the  aerodynamic  flow  parameters  in' 
the  test  section  of  a  shock  tunnel,  it  is  necessary  to  determine  the 
duration  of  usable  test  flow  as  limited  by  the  arrival  of  the  test  gas 
(driven  tube  gas)  -  driver  gas  interface.  This  investigation  will  deal 
with  a  helium-air  interface. 

A  diagram  of  the  shock  tunnel  is  shown  in  Figure  1.  A  high 
velocity,  low  density  test  flow  is  generated  by  precharging  the  driven 
tube  to  a  pressure  which  will  yield  prescribed  test  section  conditions. 
The  driver  section  is  then  pressurized  until  a  diaphragm  separating  the 
two  sections  is  ruptured  which  initiates  a  compression  shock  in  the 
driven  tube  gas.  When  this  pressure  rises  sufficiently,  a  second 
diaphragm  between  the  driven  tube  and  the  previously  evacuated  dump  tank 
ruptures,  causing  a  high  velocity,  low  density  flow  to  pass  through  the 
test  section  area.  The  parameters  of  pressure,  diaphragm  thickness, 
tube  lengths,  and  nozzle  dimensions  are  predicted  from  shock  tunnel 
practice.  This  thesis  deals  only  with  the  analysis,  development,  and 
evaluation  of  the  electron  beam  interface  detection  system.  For  a  dis¬ 
cussion  of  the  shock  tunnel  parameters,  see  Hertzberg  and  others  (1).* 

^■Numbers  in  parentheses  refer  to  similarly  numbered  references 
in  the  bibliography. 
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Fig.  1  Diagram  of  a  Shock  Tunnel 
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There  are  at  least  three  requirements  which  any  interface  de¬ 
tection  technique  must  satisfy.  First,  the  predicted  usable  test  flow 
time  is  on  the  order  of  a  few  hundred  microseconds;  therefore,  the 
method  selected  must  have  the  capability  of  tens  of  microseconds 
response  time.  Second,  it  is  not  known  whether  the  driver  gas  arrives 
as  a  well  defined  surface  or  if  there  is  a  region  of  diffused  mixing 
comparable  to  the  test  time;  therefore,  the  interface  detection  system 
must  indicate  the  time  resolved  percentage  concentration  of  driver  gas. 
Finally,  it  is  desirable  to  be  able  to  measure  the  flow  concentration 
at  a  given  position  in  the  flow;  therefore,  the  technique  should  have 
the  capability  of  nearly  point  resolution. 

Kaegi  and  Muntz  (2)  discuss  the  above  points  for  the  conditions 
of  their  study  of  interface  arrival.  Their  work  consisted  of  monitor¬ 
ing  the  density  of  nitrogen  and  helium  as  measured  by  the  electron  beam 
technique  described  by  Muntz  and  Marsden  [3) =  The  interface  was  iden¬ 
tified  with  a  decrease  in  nitrogen  density  or  an  increase  in  helium 
density.  The  two  gases  were  monitored  independently  on  separate  tunnel 
runs. 

Copper,  Miller,  and  Hameetman  (4)  and  Bird,  Martin,  and  Bell  (5) 
obtained  data  on  both  helium  and  hydrogen  driver  gas  interface  arrival. 
Their  method  consisted  of  four  fast  acting  sampling  probes  which  were 
open  at  four  different  periods  over  about  ten  milliseconds  of  apparent 
flow  time.  The  total  open  time  was  about  two  milliseconds.  In  this 
manner  four  samples  representing  average  mixtures  were  obtained.  These 
samples  were  subsequently  analyzed  by  gas  chromatographic  methods. 
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This  work  was  undertaken  because  a  valid  measure  of  run  time  is 
needed  and  the  earlier  data  were  not  considered  sufficiently  accurate 
or  complete  for  making  predictions  of  shock  tunnel  performance. 

Chapter  II  gives  an  introduction  to  electron  beam  diagnostics 
and  a  description  of  the  electron  beam  diagnostic  system  used  to  take 
the  data  in  Chapter  IV.  A  theoretical  expression  for  the  relative 
intensities  of  electron  beam,  excited  radiation  in  helium-air  gas  mix¬ 
tures  is  derived  in  Chapter  III.  Experimental  data  for  both  static 
and  transient  helium-air  gas  mixtures  are  presented  in  Chapter  IV. 
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CHAPTER  li 

ELECTRON  BEAM  GAS  DIAGNOSTICS 


I.  THE  ELECTRON  BEAM  GAS  DIAGNOSTIC  TECHNIQUE! 


A  totally  satisfactory  method  for  determining  the  unperturbed 
parameters  in  a  high  velocity  (30,000  feet/sec),  low  density  (approxi¬ 
mately  five  Torr  at  300°K  maximum)  aerodynamic  flow  has  not  been 
offered.  Although  tedious,  spectroscopic  techniques  are  particularly 
attractive  since  they  can  be  employed  without  causing  any  disturbance 
in  the  gas  during  the  time  of  measurement.  Both  emission  and  absorp¬ 
tion  spectroscopy  can  be  considered.  However,  at  the  conditions  given 
above  the  gas  does  not  emit  or  absorb  a  sufficient  amount  of  radiation 
for  quantitative  measurements.  , 

The  high  energy  (10-50  KV)  electron  electron  beam  holds  some 
promise.  By  directi.i?  a  well-collimated  beam  of  electrons  into  the 
flow  the  gas  is  caused  to  fluoresce  within  the  beam  volume.  Inelastic 
collisions  of  the  electrons  with  the  gas  molecules  result  in  the 
excitation  of  the  gas  molecules.'  In  many  common  gases  such  as  nitrogen, 
helium,  argon,  carbon  dioxide,  hydrogen,  and  others  some  of  the  energy 
resulting  from  subsequent  spontaneous  de-excitations  of  the  molecules 
is  radiated  in  the  easily-handled  optical  spectrum  [approximately  2000 

o  o  .  _  o 

A  (one  Angstrom  A  equals  10-10  meters)  to  104  A  wavelength].  Spectro¬ 
scopic  analysis  of  this  radiation  has  been  used  to  make  measurements 
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of  several  flow  parameters. 

Schumaker  and  Gadamer  (6)  used  the  electron  beam  to  measure 
nitrogen  density  below  100  millitorr  by  detecting  the  change  in  beam 
radiation  intensity  with  density  at  a  known  beam  current.  They  found 
a  linear  relationship  in  this  range.  A  we 11- collimated  beam  and  an 
accurate  beam  current  measurement  are  critical  to  this  technique. 

Muntz  (7)  used  the  relative  line  intensity  method  described  by 
Hertzberg  (8)  to  measure  the  vibrational  and  rotational  temperature  of 
the  unperturbed  gas.  Since  the  relaxation  times  for  rotational  and 
translational  temperatures  are  comparable,  this  technique  offers  a 
means  of  obtaining  the  temperature  for  use  in  thermodynamic  and  aero¬ 
dynamic  theoretical  analyses.  Measurements  of  vibrational  temperature 
provide  a  determination  of  the  degree  of  non-equilibrium. 

Muntz  (9)  made  use  of  the  Doppler  broadening  of  the  beam  excited 
emission  lines  to  obtain  data  on  the  molecular  velocity  distributions 
in  a  low  density  gas.  This  was  done  by  measuring  the  Doppler  width  of 

O 

the  5016  A  helium  line. 

The  electron  beam  therefore  provides  a  means  of  exciting  a 
selected  portion  of  a  low  density,  high  velocity  flow  in  such  a  manner 
that  optical  spectroscopic  techniques  can  be  used  to  make  measurements 
of  the  undisturbed  free  stream  parameters  of  the  flow. 

The  electron  beam  flow  diagnostic  technique  has  several  pitfalls. 
The  successful  application  of  the  electron  beam  excitation  technique 
to  any  particular  flow  field  investigation  depends  on  an  approximate 
knowledge  of  the  flow  parameters  and  the  specific  operational 
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limitations.  The  electron  beam  generating  system  must  provide  a  well 
collimated,  stable  stream  of  electrons  before  conclusions  can  be  drawn 
with  regard  to  particular  regions  of  the  flow.  This  is  especially  true 
in  density  measurements.  Care  must  be  taken  that  the  electron  beam 
hardware  and  optics  do  not  produce  disturbances  near  the  point  of 
measurement.  An  often  encountered  problem  is  background  luminosity  due 
to  flow  contamination.  Valid  data  can  only  be  obtained  when  the 
radiation  accepted  by  the  spectrometer  from  sources  other  than  the 
beam  are  reduced  to  an  acceptable  level.  The  level  of  acceptable  back¬ 
ground  radiation  is  directly  dependent  on  the  absolute  intensity  of 
beam  radiation.  In  general,  the  higher  the  density  and  beam  current 

t 

the  greater  the  background  intensity  that  can  be  tolerated. 

Another  area  which  must  be  considered  is  the  degree  to  which  the 
flow  conditions  and  the  assumed  excitation  and  de-excitation  mechanisms 
due  to  the  electron  beam  comply  with  the  actual  conditions.  Muntz, 

Abel,  and  Maquire  (10),  Muntz  (7),  and  Petrie  (11)  give  a  quite  exhaust¬ 
ive  commentary  on  this  topic. 

II.  THE  ELECTRON  BEAM  SYSTEM 

The  development  of  an  electron  beam  generating  system  and 
associated  optical  detection  system  poses  a  myriad  of  problems.  A 
diagram  of  the  system  used  for  taking  data  in  the  eight  foot  diameter 
test  section  of  Tunnel  "J"  at  Arnold  Engineering  Development  Center  is 
shown  in  Figure  2.  A  television  type  electron  gun  with  magnetic  focus¬ 
ing  was  used  to  generate  the  electron  beam.  An  early  attempt  with 
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Fig.  2  View  of  Electron  Beam  S/stem  in  Test  Section,  Looking  Downstream 
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electrostatic  focusing  was  terminated  due  to  the  difficulty  of  design¬ 
ing  a  power  supply  with  the  cathode  operated  at  high  voltage  rather 
than  near  ground  as  in  conventional  circuits.  These  electron  guns  are 
capable  of  operating  at  20  KV  while  delivering  four  to  five  milliam- 
peres  of  cathode  current.  The  filaments  are  designed  for  six  volt 
excitation  but  they  were  routinely  operated  at  eighteen  volts.  While 
this  reduced  the  lifetime  of  the  gun  it  provided  additional  emission 
and  allowed  operation  at  pressures  above  the  normal  operating  range. 

In  order  to  provide  a  steady  beam  care  must  be  taken  to  provide  low 
ripple  direct  current  power  supplies  throughout  the  electron  beam 
generation  system. 

It  is  necessary  that  the  pressure  in  the  electron  gun  be  3  x 
10-4  Torr  or  less  for  acceptable  operation.  In  order  to  maintain  this 
pressure  while  the  pressure  in  the  test  section  may  be  as  high  as  one 
Torr  for  a  few  milliseconds,  it  is  necessary  to  use  a  dynamic  pumping 
system.  This  technique  requires  firing  the  beam  through  a  small  orifice 
across  which  the  necessary  pressure  differential  can  be  maintained  by 
continuous  pumping.  The  system  in  Figure  2  uses  a  0.050  inch  diameter 
orifice  with  a  sixty  liter  per  second  diffusion  pump.  A  pressure  of 
150  millitorr  can  be  maintained  easily  in  a  small  test  section  enclos¬ 
ing  the  electron  gun  exit  orifice  with  the  electron  beam  in  operation. 
Satisfactory  operation  has  been  obtained  with  transient  flows  of  a  few 
milliseconds  in  duration  up  to  densities  equivalent  to  700  microns  at 
room  temperature. 

Passing  a  major  portion  of  an  electron  beam  which  is  one  to  two 
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millimeters  in  diameter  through  a  0.050  inch  (1.25  mm)  diameter  orifice 
places  some  rather  stringent  requirements  on  beam  alignment.  Static 
alignment  of  the  beam  is  accomplished  by  varying  the  currents  in  an  x-y 
deflection  coil  and  maximizing  the  collected  current.  Dynamic  realign¬ 
ment  is  done  by  a  feedback  control  device  developed  by  Froedge  (12) 
which  detects  current  resulting  from  misalignment  of  the  beam  and  uses' 
a  voltage  proportional  to  the  time  integral  of  this  current  to  change 
the  proper  deflection  coil  current  in  the  direction  necessary  to  re¬ 
align  the  beam.  This  control  system  has  calculated  response  time  of 
fifty  microseconds. 

Radiation  from  the  beam  is  intercepted  by  a  f/1.5  lens  and  trans¬ 
mitted  outside  the  test  section  by  a  series  of  collimating  lenses.  At 
the  partial  beam  splitter  the  light  is  separated  into  two  beams.  One 

O 

beam  is  directed  through  a  fifty  A  half  width  interference  filter 

e 

centered  on  the  5016  A  atomic  spectral  line,  which  is  prominent  in 

helium  electron  beam  excited  radiation  [see  Sebacher  (13)].  The  second 

0 

beam  is  directed  through  a  seventeen  A  half  width  interference  filter 

O 

centered  at  the  4278  A  bandhead  of  the  (0,1)  vibrational  band  in  the 
first  negative  system  of  N*,  which  is  prominent  in  nitrogen  electron 
beam  radiation  [see  Sebacher  (13)].  Both  radiation  outputs  are  detected 
with  high  gain  RCA  7650  photomultipliers  whose  outputs  are  in  turn 
displayed  on  oscilloscopes.  In  addition  the  photomultiplier  outputs  are 
fed  into  a  diode  type  quarter-square  multiplier  connected  as  a  divider. 
This  device  is  discussed  in  Chapter  IV.  In  this  way  the  "in  line"  ratio 
of  helium  intensity  divided  by  air  (nitrogen)  intensity  is  obtained. 
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A  tungsten  strip  lamp  was  used  as  a  radiation  intensity  stand¬ 
ard.  .  Radiation  from  the  tungsten  strip  lamp  was  directed  into  the 
detector  portion  of  the  optics  by  a  clear  glass  beam  splitter.  While 
the  purpose  of  the  standard  lamp  was  to  monitor  the  gain  factors  of 
the  optical  system  over  the  several  weeks  that  data  were  taken  the 
calibration  can  also  be  used  to  make  estimates  of  absolute  intensity 
sensitivity. 

A  flat  catcher  plate  about  twelve  inches  in  diameter  and  eight¬ 
een  inches  from  the  beam  exit  orifice  was  used  to  collect  the  electron 
beam.  While  it  definitely  did  not  provide  a  suitable  measure  of  beam 
current  for  density  measurements,  it  did  give  a  measure  of  beam  quality 
and  continuity  during  the  transient  flow. 
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CHAPTER  III 

CALCULATION  OF  EXPECTED  INTENSITIES 


I.  ESTIMATE  OF  PHOTON  YIELD  FROM  ELECTRON  BEAM  INTERACTION  WITH  A  GAS 


The  general  equation  for  the  intensity  of  a  spectral  line,  in 
emission  is  given  in  Herzberg  (8): 


I--  ■  N  h  v__  A__. 
nm  n  ran  iim 


CD 


In  this  equation  1^  is  the  intensity  of  the  radiation  emitted  into  4tt 
steradians  when  the  molecules  or  atoms  of  a  gaseous  ensemble  decay  from 
energy  state  n  to  energy  state  m,  Nn  is  the  number  of  molecules  or  atoms 
in  energy  state  n,  is  the  frequency  of  the  light,  h  is  Planck's 
constant,  and  A^,  is  the  Einstein  probability  of  a  spontaneous  emission. 
Values  of  A^,  are  available  for  most  of  the  prominent  emission  lines  in 
the  common  gases.  The  frequency  of  the  light  can  be  determined  direc¬ 
tly.  It  remains  to  determine  the  number  of  molecules  or  atoms  in  the 
excited  state  n  as  a  result  of  electron  beam  excitation. 

Consider  the  energy  level  model  in  Figure  3.  The  following  con¬ 
ditions  will  be  assumed  throughout  the  remainder  of  this  analysis: 

1.  Only  one  lower  or  ground  level  l  will  be  considered.  This 
is  essentially  true  up  to  a  few  hundred  degrees  Kelvin  above 
room  temperature  for  most  gases  since  the  population  of 
higher  states  is  very  low. 
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2.  The  electron  beam  excitation  is  a  direct  excitation.  There 
is  no  cascade  population  of  the  n  excited  energy  level  by 
higher  levels.  Also  there  are  no  excitations  with  inter¬ 
mediate  stops  at  lower  energy  levels. 

3.  There  are  no  collisional  de-excitations  from  the  excited 
state  n. 

The  equation  for  the  number  of  molecules  or  atoms  being  excited 
to  the  state  n  per  unit  volume  and  per  unit  time  may  be  written. 


N 


ex 


e  Av 


(2) 


molecules 

excited 


sec  •  nr 


J  coul/sec 
coul 

. 


m 


<  \ 

) 

H 

2 

m 

>< 

molecules  l 

3 

m 

k  j 

‘  i 

where  is  the  beam  current,  e  is  the  electronic  charge,  Aj,  is  the 
area  of  the  beam,  Q£n  is  the  cross-section  for  excitation  from  state 
£  to  state  n,  and  is  the  number  of  molecules  per  unit  volume  in 
state  £.  The  equation  for  the  number  of  molecules  or  atoms  de-excited 
from  state  n  may  be  written. 


ndx  *  T.  W  ts> 

m 

o 

The  decay  times  for  spontaneous  emission  are  on  the  order  of  10 
seconds  for  electric  dipole  transitions.  Within  a  few  decay  times 
after  the  beam  is  turned  on  the  number  of  excitations  per  second  must 
equal  the  number  of  de-excitations  per  second. 
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e  Av 


Q 


Jin 


ni  =  Y. 

m 


V  Nn  • 


(4) 


In  this  steady  state  condition  Nn  is  a  constant.  Also  ^  Anm  =  l/tn, 

m 

where  in  is  the  mean  lifetime  of  state  n  as  defined  in  Herzberg  (8) . 


Nn  = 


!b  Q&n  NJl  Tn 
e  A u 


(5) 


Substituting  this  result  into  Equation  1  gives 


t  -  h  vnm  Aim  *b  Q&n  Nft  Tn  (6) 

nm  “  e  Ab 

The  above  expression  gives  the  energy  per  unit  volume  per  second 
radiated  isotropically  into  4tt  steradians  due  to  the  electron  beam 
excitation  of  a  single  energy  transition  or  a  group  of  transitions 
which  can  conveniently  be  considered  as  a  single  transition.  It  is  of 
interest  to  note  that  if  it  were  valid  to  consider  only  one  de-excita¬ 
tion  transition  the  product  of  A^  and  xn  would  be  one.  In  that  case 
these  two  parameters  would  not  appear  in  the  expression  for  intensity. 

In  order  to  determine  the  absolute  intensity  available  to  a 
photodetector  it  is  necessary  to  consider  the  geometry  of  the  situation. 
Assuming  an  isotropic,  homogeneous,  optically  thin  medium  the  total 
power  If  radiated  into  4tt  radians  will  be  the  product  of  the  power 
radiated  per  unit  volume  and  the  volume  of  beam  excited  gas  observed. 

If  =  *nm  ha  ^ 

where  Jl^  is  the  length  of  beam  observed.  Substituting  Equation  6  into 

Equation  7, 
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I  b  h  ^  ^nm  Ib  Qtn  Ni*b  Tn  ,  rgj 

T 

The  steradiancy  intercepted  by  an  optical  system  can  be  related  to  the 
system  f  number. 


f  = 


Where  f  is  the  system  f  number,  D  is  the  lens  stop  diameter,  and  is 
the  focal  length  of  the  optical  component  in  the  system  with  the  maxi¬ 
mum  f  number.  The  intercepted  solid  angle  ft  is 


ft  =  A_  .  (10) 

d2 

Where  A  is  the  collector  lens  area,  and  d  is  the  distance  from  the 
source.  Clearly,  if  d  is  made  equal  to  the  focal  length. 


ft  = 


(ID 


The  total  amount  of  radiation  intercepted  I^n  can  be  determined  from 
Equations  8  and  11: 


r.  „  k  vnm  Anm  Ib  *5in  *b  Tn 

iln  - = - 

16  fz  e 


(12) 


Introducing  a  factor  K  to  allow  for  the  optical  system  losses  in  len¬ 
ses,  mirrors,  and  filters  gives  the  intensity  IpD  incident  on  a 
photodetector. 

K  h  Vnjj,  Aj^,  Ijj  Q4n  rn 


PD 


16  f2  e 


(13) 


In  general,  the  emission  intensity  is  a  function  of  beam  current 
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and  gas  density.  The  experiments  of  Schumacher  and  Gadamer  (6)  and 
Petrie  (12)  show  that  the  intensity  is  linear  with  beam  current  and  gas 
density  up  to  at  least  one  milliampere  and  a  few  hundred  millitorr  of 
pressure  at  room  temperature  in  air.  Therefore,  within  the  above 
limitations  the  intensity  can  be  considered  linear  with  density. 

Introducing  the  ideal  gas  law  in  the  form  shown  by  Kieffer  and 
Dunn  (14),  parameters  more  easily  measured  experimentally  are  obtained. 


JPD  = 


9.62  x  10*4  K  h  Vjm  Ij,  Q^n  P  *{,  'n 


(14) 


16  f  e  T 

where  P  is  the  pressure  in  Torr,  T  is  the  temperature  in  °K,  and  the 
other  parameters  are  consistent  with  mks  units. 

The  most  inaccurate  parameters  in  Equation  14  are  the  excitation 
cross-section  Qjn,  the  transition  probability  A^,  and  the  lifetime  of 
the  excited  state  tn.  The  parameters  and  tn  have  been  measured 
within  an  order  of  magnitude  while  Qtn  may  be  more  than  an  order  of 
magnitude  in  error.  Therefore,  the  accuracy  of  the  calculation  under 
the  assumptions  rests  primarily  on  the  accuracy  of  the  excitation 
cross-section  Q4n. 

A  theoretical  value  for  the  excitation  cross-section  can  be 
obtained  by  applying  the  Bom  approximation  for  inelastic  excitation 
of  molecules  by  fast  electrons.  For  incident  electron  energies  much 
in  excess  of  the  excitation  threshold  of  the  n  level  the  collision 
cross-section  for  optically  allowed  excitation  transitions  to  the  n 
state  from  the  ground  state  may  be  expressed  by  the  following  equation 
from  Mott  and  Massey  (15). 
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.  64tt^ 

(2t  n  v/h)2h4 


(15) 


Where  jZ^!  is  the  matrix  element  of  the  transition,  En  is  the  energy 
of  the  excited  state,  is  the  energy  of  the  ground  state,  v  is  the 
velocity  of  the  incident  electron,  m  is  the  mass  of  the  electron,  h  is 
Plank's  constant,  and  e  is  the  electronic  charge.  Prom  this  express¬ 
ion  it  is  clear  that  the  cross-section  for  an  optically  allowed  excit¬ 
ation  transition  varies  as  Ee”*  in  Eg,  where  Eg  is  the  energy  of  the 
incident  electron-.  The  functional  dependence  of  the  excitation  cross- 
section  on  electron  energy  shown  in  Equation  15  will  be  used  later  to 
estimate  the  cross-section  for  10  Kev  electrons  from  experimental  data 
for  low  energy  electrons., 


II.  INTENSITY  RELATIONS  IN  HELIUM-AIR  MIXTURES 


It  is  necessary  to  consider  specifically  the  radiation  resulting 
from  the  interaction  of  the  beam  with  helium  and  air.  The  energy  level 
diagrams  and  prominent  transitions  for  helium  and  nitrogen  due  to  high 
energy  electron  beam  excitation  are  given  in  Figure  4. 

When  a  10-50  KV  electron  beam  passes  through  helium  at  densities 
below  a  few  Torr  and  room  temperature,  the  most  prominent  radiation  is 

9 

the  atomic  line  at  5016  A.  This  radiation  results  from  excitation  of 
the  l'S  ground  state  of  the  helium  atom  to  the  3'P  excited  state  by  the 

O 

fast  electrons  and  the  subsequent  spontaneous  de-excitation  of  5016  A 
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energy  to  the  2'S  excited  state.  In  the  following  discussion  the 

O 

parameters  necessary  for  calculating  the  expected  intensity  of  5016  A 
radiation  will  be  obtained  in  the  light  of  the  assumptions  made  in 
deriving  the  equation  for  intensity  of  radiation  from  electron  beam 
interaction  with  a  gas. 

Pirst,  consider  the  validity  of  assuming  a  single  ground  state. 
The  ratio  of  the  number  of  molecules  in  excited  states  to  the  number 
in  the  ground  state  is  given  by  the  Boltzman  factor  as  shown  by 
Herzberg  (16). 

N„  gn  -(En-Et)/kI 

"7-*7e  •  tl 

where  gft  and  g4  are  the  statistical  weights  of  the  upper  and  lower 
states,  respectively,  and  the  other  parameters  are  as  defined  pre¬ 
viously.  For  the  l'S  ground  state  and  2'S  excited  state,  which  is  the 
lowest  excited  state  of  helium,  Equation  16  becomes. 


N  2.35  x  105 

/  *  e  T  .  (17) 

Equation  17  shows  that  for  temperatures  below  50, 000° K,  99  per  cent  of 
the  molecules  are  in  the  l'S  ground  state.  Also  the  probability  of 
excitation,  which  is  comparable  to  the  transition  probability  in  a 
relative  sense,  from  excited  states  to  the  3'P  level  is  less  than  that 
from  the  l'S  ground  state.  Therefore,  only  one  lower  level  need  be 
considered  and  the  single  level  assumption  is  valid. 

In  order  to  calculate  the  excitation  cross-section  for  10  KV 
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electrons,  the  value  from  Bauer  and  Bartky  (17)  for  the  excitation  of 
the  3'P  level  by  150  volt  electrons,  which  is  well  above  the  excitation 
level  where  the  Bom  approximation  and  experimental  values  agree,  will 
be  used  in  conjunction  with  the  Ee’*  in  Ee  electron  energy  dependence. 


*£n_ 

<*150 


Ee"l  tn  Ee 
E150  ln  E1S0 


(18) 


The  excitation  cross-section  for  10  KV  electrons  calculated  from 
Equation  18  is  0.82  x  lO”2^  m2. 

The  population  of  the  3'P  state  is  due  to  direct  excitation  since 
transitions  to  higher  states  (from  which  cascading  could  occur)  and 
lower  states  (an  intermediate  step  in  two  step  excitation)  are  not 
optically  allowed.  This  can  be  easily  deduced  from  the  table  of  tran¬ 
sition  probabilities  given  by  Wiese,  Smith,  and  Glennon  (18). 

There  are  two  major  factors  to  consider  for  de-excitation  of 
the  3'P  state.  For  a  first  order  calculation  collisional  de-excitation 
effects  will  not  be  considered.  However,  at  high  density  this  would 
be  a  major  contributor.  The  other  factor  is  resonant  excitation  - 
de-excitation  of  the  3’P  state. 

The  3'P  excited  level  has  a  large  cross-section  for  de-excitation 


to  the  l'S  ground  state  emitting  a  537  A  photon.  This  537  A  photon  is 
quite  readily  absorbed  by  unexcited  helium  atoms,  and  can  escape  from 
an  enclosure  only  after  many  absorptions  and  re-emissions.  Clearly 
this  phenomena  will  have  an  effect  on  the  lifetime  and  the  population 
of  the  3'P  level  of  helium.  As  previously  stated  any  diagnostic 
technique  of  use  in  shock  tunnels  must  have  the  capability  of  a 
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response  time  of  a  few  tens  of  microseconds.  It  is  necessary  then  to 
investigate  the  effect  of  resonance  diffusion  on  the  decay  time  of  the 

e 

5016  A  helium  line.  The  following  analysis  of  the  effect  of  resonance 
radiation  on  excited  level  lifetimes  will  follow  the  work  of  Holstein 
(19). 

If  N0  atoms  are  excited  to  a  reference  level  prior  to  t  =  0,  then 
the  number  remaining  in  this  level  as  a  function  of  time  is  determined 
by  the  equation, 

N  =  N0  e  ”qAt  ,  (19) 

where  A  is  the  probability  per  unit  time  for  emission  (the  Einstein 
Coefficient)  of  a  quantum  of  resonance  radiation  from  an  isolated 
atom,  and  q  is  the  imprisonment  factor.  The  imprisonment  factor  is  less 
than  unity  and  depends  on  the  density,  the  shape  of  the  containing 
vessel,  and  the  dimensions  of  the  containing  vessel.  The  magnitude  of 
q  varies  from  zero  to  one.  When  the  gas  density  is  low  q  is  equal  to 
one  and  the  lifetime  of  the  3'P  level  is  determined  entirely  by 
spontaneous  transition  probabilities.  At  higher  gas  densities  q  goes 
to  zero  and  the  effect  of  resonance  trapping  is  significant. 

Heron,  McWhorter,  and  Roderick  (20),  using  pure  helium,  have 
shown  that  the  lifetime  of  the  3'P  level  is  determined  largely  by  the 
phenomenon  of  imprisonment  of  resonance  radiation,  is  pressure- 
dependent,  and  is  of  the  same  order  of  magnitude  as  the  lifetimes 

O  O 

measured  for  the  3914  A  and  4278  A  ionized  nitrogen  bands  in  a 
mixture  of  helium  and  nitrogen.  J 

i 

The  Einstein  Coefficient  A  is  the  reciprocal  of  the  mean  life¬ 
time  of  the  resonance  level  in  the  absence  of  imprisonment.  The 
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imprisonment  factor  q,  as  given  by  Holstein  (19)  for  a  cylinder  of 
infinite  length  with  the  radiation  source  at  the  centerline,  is 


1.60 _ 

k0R  (it  Zn  k0R)l/2 


(20) 


where  R  is  the  radius  of  the  containing  vessel  and,  for  a  Doppler- 
broadened  line,  kQ  is 


k  -  ill 

o  "  8  ir 


£2 

*1 


*n  Vo  * 


1/2 


(21) 


Where 


is  the  most  probable  thermal  speed  of  a  helium  atom,  \  is  the  wave¬ 
length  of  the  radiation,  N  is  the  number  of  unexcited  helium  atoms  per 
cubic  centimeter,  k  is  Boltzman's  constant,  m  is  the  molecular  mass, 
and  gj  and  g2  are  the  statistical  weights  of  the  resonance  and  ground 
levels,  respectively. 

As  previously  stated  any  diagnostic  technique  of  use  in  shock 
tunnels  must  have  the  capability  of  a  response  time  of  a  few  tens  of 
microseconds.  It  is  necessary  then  to  investigate  the  effect  of 

O 

resonance  diffusion  on  the  decay  time  of  the  5016  A  helium  line. 

In  the  case  of  the  electron  beam,  the  excited  states  are  pro¬ 
duced  at  a  rate  Q0  per  unit  time,  the  rate  depending  mainly  on  the 
electron  current.  The  total  excited  in  dt  is  Q0dt.  The  intensity  of 

O 

the  5016  A  line  can  be  expressed  as  a  function  of  time  as 
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I5016  ^  m  h  VS016  N3’P  W  A3'P-2'S  * 


(23) 


The  incremental  number  of  atoms  remaining  in  the  resonant  state  due  to 
the  number  of  atoms  excited  in  interval  dt  is 

-q'A't  -q'A't 

dN  =  Q0  dt  e  =  Q0  e  dt.  (24) 


Integrating 

N(t)  = 


(25) 


0  -q'A’t 

«W  “  Cl  -  •  )•  (26) 

Micro  q'A'  is  the  sum  of  all  transition  probabilities  and  products  of 
imprisonment  factors  and  transition  probabilities. 

In  the  case  of  the  3'P  state  q'A'  involves  primarily  the  5016 

o  o  o 

A  transition  and  the  537  A  transition,  but  only  the  537  A  transition  is 
resonant  and  has  an  imprisonment  factor;  therefore. 


q’A'  =  q  A537  +  A5016  =  _L_  .  (27) 

l3'P 

Thus, 


N31P  (t)  = 


q  A537  +  A5016 


Substituting  Equation  28  into  Equation  23, 


'it»W  -  h  "5016  q-^5°16 

q  a537 


016  . 

+  A  1 

5016 


e  " (qA537+A5016^ j _  (28) 
e  'CqAS37*AS016>tl 
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For  example,  assume  the  following  representative  values  taken  from 
Heron,  McWhorter,  and  Roderick  (20). 


'5016  s  **55  x  sec"*. 

(30) 

o  _  1 

.537  =  5.6  x  10  sec  , 

(31) 

5016  =  •"  x  10*5  Hertz. 

(32) 

Substituting  into  Equation  29, 

Q 

I5016(t)  =  56q  +  1.35 

As  pressure  goes  from  zero  to  infinity,  q  goes  from  unity  to 
zero.  The  extreme  effect  of  resonance  trapping  occurs  when  the  pres¬ 
sure  approaches  infinity.  The  decay  time  in  this  limiting  case  is  the 

O 

decay  time  of  the  5016  A  transition  which  is  0.074  microseconds. 
Therefore,  resonance  diffusion  will  not  preclude  use  of  this  line  for 
data  taken  in  observation  times  greater  than  one  microsecond. 

O 

Finally,  considerable  survey  spectra  have  shown  the  5016  A  line 
sufficiently  separated  from  other  radiation  in  helium  and  air  to  be 
quite  acceptable  for  diagnostic  applications. 

The  above  parameters  and  considerations  compiled,  the  expected 
intensity  according  to  Equation  14  may  be  calculated.  Repeating  for 
clarity, 

V3'P-2'S  =  =  *  1014  sec  *,  (34) 

where  c  is  the  speed  of  light, 


(56q  +  1.35)xl07t 


(33) 


25 


AEDC-TR-68  -88 


^3 iP-2 * 5  =  1.35  x  10  sec 

(35) 

QltS_3.p  =  0.82  x  10'23  m2. 

(36) 

q  =  1, 

(37) 

and 

T  =  1  =  1.75  x  10“9  sec. 

nmin  "537  +  A5016 

(38) 

Substituting  into  Equation  14, 

Ij^  «  1.15  x  10 


.10  9.62  x  1024  KHeh  IfeP  *b 


watts. 


16  fz  e  T 


where  P  is  in  Torr,  T  is  in  °K,  is  the  transmission  factor  for  the 
optical  system,  and  the  other  parameters  are  consistent  with  mks  units. 


In  a  similar  manner  the  interaction  of  an  electron  beam  with 


air  will  be  discussed.  In  air,  under  conditions  comparable  to  those 
considered  for  helium,  the  most  prominent  radiation  is  due  to  what  is 


called  the  first  negative  system  of  ionized  nitrogen.  This  system  con¬ 
sists  of  several  vibrational  bands  with  the  strongest  band  having  its 
bandhead  at  3914  A.  This  system  results  from  excitation  of  the 
ground  electronic  state  of  the  stable  nitrogen  molecule  to  the  B2£ 
excited  state  of  the  singly  ionized  nitrogen  molecule  by  the  fast 


electrons  and  subsequent  spontaneous  de-excitation,  which  results  in 
the  first  negative  system,  to  the  X^  Y.  state  of  the  ionized  mole¬ 
cule.  The  nitrogen  system  is  inherently  much  more  complicated  than 

O 

atomic  helium  since  each  vibrational  band  is  approximately  forty  A 
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wide  and  consists  of  a  number  of  rotational  lines  [see  Muntz  (7)].  How 
ever,  by  considering  each  vibrational  band  to  be  a  single  line,  as  it 
would  appear  on  a  spectrograph  with  poor  resolution,  an  analysis 
similar  to  the  analysis  of  helium  may  be  carried  through. 

Consider  the  validity  of  assuming  a  single  ground  vibrational 
state.  The  ratio  of  the  number  of  molecules  in  the  first  excited 
vibrational  state  of  the  neutral  molecule  to  the  number  in  the  ground 
state  is  given  by  Herzberg  (8)  as  1.4  x  10“^  at  300°K  and  0.035  at 
1000°K.  At  a  temperature  of  1000*K,  96.5  per  cent  of  the  molecules  are 
in  the  ground  vibrational  state.  This  requirement  is  not  critical 
since  the  excitation  cross-section  used  for  subsequent  calculations 

e 

is  the  total  cross-section  for  the  excitation  of  the  4278  A  vibrational 
band. 

A  rigorous  calculation  of  the  excitation  cross-section  using 
the  Born  approximation  has  been  done  by  Homkohl  (21).  These  cal¬ 
culations  have  been  done  for  the  purpose  of  studying  high  energy 
electron  excitation  and  have  shown  good  functional  agreement  with 
experimental  data  by  Sheridan  (22).  These  data  will  be  used  for 
nitrogen  excitation  cross-section. 

A  considerable  amount  of  investigation  on  the  excitation- 
emission  path  in  nitrogen  has  been  done  by  several  authors.  Muntz  (7) 
gives  an  exhaustive  discussion  with  the  conclusion  that  the  excita¬ 
tion  to  the  state  is  a  direct  excitation  and  it  may  be  con¬ 

sidered  the  only  path  of  excitation.  Measurements  of  the  mean  life¬ 
time  by  Bennet  and  Dalby  (23)  have  shown  that  it  is  not  affected  by 
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density  up  to  at  least  1.5  x  1016  atoms  per  cm3  (470  millitorr  at  room 

temperature).  Therefore,  this  system  should  be  satisfactory  for  air 

diagnostics  and  the  model  assumed  previously  should  be  acceptable  for 

first  order  absolute  intensity  calculations.  Finally,  as  was  done  with 

helium,  considerable  search  spectra  were  taken  to  determine  what  por- 

tion  of  the  system  should  be  used.  Although  the  3914  A  band  is  the 

0 

strongest  the  4278  A  band  was  chosen  because  early  experiments  demon* 

O 

strated  a  noticeable  interference  from  the  3889  A  atomic  line  in  helium. 
The  above  parameters  and  considerations  compiled,  the  expected 

O 

intensity  according  to  Equation  14  may  be  calculated  for  the  4278  A 
(0,1)  band  of  nitrogen.  The  parameters  are, 


and 


-  Njx2! 

(0,1)  - 

c/4278  A  =  7.0  x  1014  sec"1. 

(40) 

-n|x2£ 

O 

%• 

h-» 

n 

0.220  x  107  sec"1. 

(41) 

Ve 

-N2B2  E 

»  0.149 

x  10"22  m2. 

(42) 

tN+b2£  f0*1) 


6.58  x  10"8  sec. 


(43) 


Substituting  into  Equation  14, 

t  i  c,  . __9  9.62  x  1024  KN2  h  Ib  PN2  *b 

IN  =  1.51  x  10  y  - - - - - = - 2 - watts.  (44) 

1  16  fz  e  T 

Where  P  is  in  Torr,  T  is  in  °K,  %2  is  the  transmission  factor  for  the 
optical  system,  and  the  other  parameters  are  in  mks  units. 

In  order  to  get  a  more  positive  indication  of  the  arrival  of  a 
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given  percentage  of  helium,  say  ten  per  cent,  at  the  test  section 

% 

where  the  flow  density  is  often  changing  quite  rapidly,  the  ratio  of 
helium  intensity  to  air  intensity  was  measured  directly.  Density  has 
a  compound  effect  on  intensity  because  intensity  is  directly  dependent 
on  density  and  beam  current  while  beam  current  is  dependent  on  density 
due  to  beam  attenuation.  The  constants  KHe  and  Kn2  in  Equations  39 
and  44  are  the  constants  for  a  particular  channel  as  shown  in  Figure  2, 
page  8. 

According  to  Dalton's  law  of  partial  pressures,  each  gas  in  a 
mixture  of  gases  exerts  pressure  in  proportional  part  by  volume. 
Therefore,  in  a  volume  percentage  mixture  of  helium  and  nitrogen  the 
pressure  in  Equations  39  and  44  become  the  partial  pressures  of  helium 
and  nitrogen.  Taking  the  ratio  of  Equation  30  to  Equation  44, 

THe  =  1-15  x  10"10  KHe  PHe 
%  =  1*51  x  lO'9  KN2  PNz 

The  parameters  of  beam  current,  beam  area,  length  of  beam  observed, 
temperature,  and  system  f  number  do  not  appear  in  Equation  45. 

Equation  45  may  be  written  in  terms  of  a  volume  percentage: 


•  0.076 


KHe  (%He) 

V5? 


(46) 


Since  nitrogen  is  seventy-eight  per  cent  by  volume  of  air,  Equation  46 
may  be  written  for  air: 


KHe  (%He) 

°-098  Ka  (%  A)  ’ 


(47) 
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The  constant  Kjje/K^  is  a  characteristic  of  the  electro-optical  system 
and  may  be  determined  using  a  standard  intensity  source.  In  practice, 
it  has  proven  convenient  to  plot  the  intensity  ratio  versus  the  per¬ 
centage  ratio  giving  a  linear  plot  with  slope  0.098  K^/K^.  In  this 
manner  the  system  calibration  curve  may  be  obtained  from  a  single 
calibration  point. 

As  a  cursory  check  on  Equation  46  the  data  of  Sebacher  (13) 
were  considered.  These  data  consist  of  spectra  taken  in  static  mix¬ 
tures  of  helium  and  air  excited  by  a  high  energy  electron  beam. 
Assuming  the  data  are  corrected  for  system  constants.  Which  appear 
to  be  the  case  since  the  ratio  of  (0,0)  intensity  to  ^(0,1) 

intensity  is  3.5  compared  to  a  theoretical  value  of  3.35,  the  ratio 
e  • 

of  5016  A  helium  intensity  to  4278  A  nitrogen  peak  intensity  is  0.96 
for  a  mixture  of  87  per  cent  helium  and  13  per  cent  nitrogen.  From 
Equation  46  a  ratio  of  0.51  is  calculated  for  the  same  mixture  and  a 
KHe/KA  of  unity.  This  agreement  is  reasonable  considering  the  order 
of  magnitude  accuracy  of  some  parameters  in  the  analysis. 
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CHAPTER  IV 

MEASUREMENTS  OF  PERCENTAGE  CONCENTRATION  IN  HELIUM-AIR  MIXTURES 


I.  RELATIVE  INTENSITY  MEASUREMENT 

As  the  previous  analysis  has  shown,  there  is  a  considerable 
advantage  in  obtaining  the  ratio  of  intensities  when  measuring  volume 
percentages  in  a  mixture  of  helium  and  air.  This  could  be  done  by 
utilizing  the  individual  photomultiplier  outputs  for  helium  and  air 
intensity.  However,  as  is  the  usual  case  for  shock  tunnel  data,  this 
information  is  recorded  on  oscilloscopes  which  makes  point  by  point 
data  reduction  inaccurate.  Also,  as  is  shown  in  Figure  5,  due  to  the 
dependence  of  absolute  intensity  on  several  variables  the  individual 
oscilloscope  traces  are  quite  irregular.  This  fact  contributes 
additional  error  to  ratio  data  obtained  from  the  individual  oscillo¬ 
scope  traces.  Therefore,  it  was  necessary  to  obtain  the  "in  line" 
ratio  as  the  data  were  generated  during  a  shock  tunnel  run. 

The  Donner  Model  3732  diode  quarter-square  multiplier  used  to 
take  the  "in  line"  intensity  ratio  is  shown  in  block  diagram  form  in 
Figure  6,  page  33.  The  heart  of  this  operation  is  a  diode  and 
resistor  network  with  an  output  proportional  to  the  square  of  the  in¬ 
put.  By  proper  selection  of  resistors  in  a  parallel  combination  of 
several  sets  of  one  resistor  and  a  diode  in  series  this  squaring 
function  can  be  closely  approximated  in  a  "piecewise  linear"  fashion. 
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Operational  Amplifier 


Fig.  6  Block  Diagram  of  Ratio  Device 
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Two  of  these  squaring  circuits  can  then  be  combined  to  solve 

h  j2-?)2  -  k2  (z-r)2  * 

Assuming  Z  in  Figure  6  to  be  the  output  of  the  operational  amplifier, 
the  multiplier  output,  according  to  Equation  48,  is  K3YZ.  Summing  the 
inputs  to  the  operational  amplifier,  realizing  the  amplifier  draws 
negligible  current,  gives 

X  *  -  K3YZ.  (49) 

Therefore, 

z  ■  -  4r  •  cso) 

If  the  helium  photomultiplier  output  is  made  the  X  input  and  the  air 
photomultiplier  output  is  made  the  Y  input  the  intensity  ratio  is 
obtained  directly. 

A  comparison  of  the  individual  photomultiplier  outputs  and  the 
ratio  output  from  the  divider  is  shown  in  Figure  5,  page  32.  The 
advantage  of  this  technique  is  evident.  Despite  the  irregularity  of  the 
individual  helium  and  air  intensities  the  ratio  was  readily  obtained. 

All  of  the  dynamic  data  were  obtained  using  the  ratio  device. 

II.  RELATIVE  INTENSITY  MEASUREMENTS  IN  A  STATIC  GAS 


In  order  to  assess  system  quality  and  determine  static  cal¬ 
ibration  constants  a  static  evaluation  was  performed.  A  small  static 
test  section  was  placed  over  the  electron  'gun  exit  orifice.  The  beam 
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was  viewed  by  the  optics  through  a  plexiglass  port.  In  this  manner  the 
system  was  calibrated  in  the  same  physical  configuration  in  which  it 
was  used  to  make  later  dynamic  tests.  Mixtures  of  helium  and  air  from 
high  pressure  bottles  were  leaked  into  the  test  section  and  the  ratio 
of'  helium  channel  output  to  aiT  channel  output  was  taken.  This  method 
of  obtaining  mixtures  was  found  to  be  much  more  consistent  than  any 
form  of  low  pressure  mixing.  The  static  calibration  data  are  plotted 
in  Figure  7.  Since  the  primary  interest  was  in  detecting  the  arrival 
of  a  ten  per  cent  concentration  of  helium  extensive  calibration  was  only 
done  up  to  mixtures  of  twenty-five  per  cent.  This  range  is  most 
accurate  for  relative  measurements  because  accurate  high  helium  per¬ 
centages  are  not  as  easily  obtained.  Leaks  in  the  vacuum  system  are 
more  negligible  at  high  percentages  of  air. 

Before  the  data  in  Figure  7  can  be  compared  with  the  theoretical 
result  in  Equation  47  it  is  necessary  to  look  more  closely  at  the 
transmission  of  the  optical  system.  A  composite  of  the  optical  system 
components  with  variable  transmission  versus  wavelength  is  shown  in 
Figure  8,  page  37.  From  Figure  8,  it  is  observed  that  the  constant 
KHe/KA  is  determined  in  the  following  manner: 


[He]  =  KHe  J  SC*>  ^  dX 
[a  'SL  Ka  f  S(A)  FA(X)  dA 


(51) 


Where  S(x)  is  the  standard  lamp  output  function  and  FHe(X)  and  Fa(a) 
aTe  the  filter  transmission  functions.  An  estimate  of  %e/KA  is 
obtained  by  simply  considering  the  standard  lamp  intensity  at  the 
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Fig.  7  Interface  Detection  System  Static  Calibration 
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filter  transmission  peak,  the  relative  peak  transmission  of  the  filters, 
and  the  relative  half  width  of  the  filters.  Using  these  approximations 
in  Equation  51, 

pHel  =  fye  ^5016  FHe(peak)  FHe (HW)  Fblue(50l6)  > 

\XA  IsL  *  kA  S4278  fA  (peak)  ^A  (HW)  Fblue(4278)  ’ 

Substituting  values  from  Figure  8, 

/xHe \  KHe  2.06  x  SO  x  50  x  3  ,c«r> 

1!A  Isl  =  KA  0,47  x  35  x-  20  x  65  1  * 

The  measured  ratio  of  IHe  to  1^  for  the  standard  lamp  was  1.2.  Sub¬ 
stituting  this  value  into  Equation  52  gives 


K^-«1.7.  (S3) 

This  interface  detection  system  was  designed  to  measure  per  cent 
concentration  of  helium.  Interpretation  of  data  with  regard  to 
absolute  determination  of  excitation  cross-sections  or  other  parameters 
is  not  reasonable.  However,  it  should  compare  on  an  order  of  magnitude 
basis  with  relative  measurements  of  helium  and  nitrogen  parameters. 

Using  the  approximation  that  the  nitrogen  filter  has  uniform  response 
throughout  the  (0,1)  band,  it  is  quite  easy  to  compare  the  prediction 
of  Equation  47  with  the  results  shown  in  Figure  7,  page  36.  Substitut¬ 
ing  the  value  for  K^g/K^  from  Equation  53  into  Equation  47  gives 


*He 


=  0.17 


%  He 
V  Air 


From  Figure  7,  page  36,  it  is  seen  that 


(54) 


IHe 

lA 


0.90 


%  He 
%  Air 


(55) 


There  are  two  primary  sources  of  error  which  contribute  to  the 
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inequality  of  Equations  54  and  55.  As  has  been  stated  the  error  in  the 
theoretical  parameters  may  easily  total  an  order  of  magnitude.  In 
addition,  the  errors  in  data  acquisition  contribute  directly  to  the 
values  in  Equation  55. 

III.  RELATING  INTENSITY  MEASUREMENTS  IN  A  TRANSIENT 
GAS  MIXTURE  OF  HELIUM  AND  AIR 


The  ultimate  purpose  of  the  interface  detection  system  is  to 
measure  the  time-resolved  percentage  of  helium  in  a  helium-air  mixture 
which  is  transient  with  respect  to  space  and  time.  Data  were  taken 
at  a  single  point  in  the  shock  tunnel  test  section. 

Two  experiments  were  undertaken  to  determine  the  validity  of 
the  system.  First,  the  driven  tube  was  charged  with  zero,  ten,  and 
twenty-five  per  cent  volume  mixtures  of  helium  and  air.  Operation  of 
the  shock  tunnel  shown  in  Figure  1,  page  2,  as  near  its  optimum  con¬ 
dition  as  possible  with  varying  mixtures  in  the  driven  tube  gave  a 
known  gas  mixture  in  the  test  section  for  from  ten  to  twenty  milli¬ 
seconds.  Actual  traces  of  the  ratio  device  output  are  shown  in  Figure 

i 

9.  The  maximum  apparent  step  response  of  two  milliseconds  indicates 
that  there  are  no  transient  effects  of  resonance  diffusion  or  collisions 
which  invalidate  the  data  shown  in  Figure  9  which  is  essentially  con¬ 
stant  for  ten  to  twenty  milliseconds.  Complete  results  of  the  dynamic 
evaluation  are  plotted  in  Figure  10,  page  41.  Each  data  point  repre¬ 
sents  a  funnel  run.  From  this  curve  the  following  dynamic  equation  is 
derived, 
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t  5  msec ./cm. 


Flow  Start  Flow  Breakdown 

Fig.  9  Oscilloscope  Traces  of  Intensity  Ratio  Versus  Time 
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Fig.  10  Interface  Detection  System  Dynamic  Evaluation 
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He  .  __  %  He 

X*  B  °-87  OIF 


(56) 


In  the  second  experiment  a  buffer  section  and  diaphragm  were 
added  so  that  the  arrival  time  of  the  helium  interface  could  be  varied 
in  a  predictable  manner.  The  results  of  this  experiment  are  shown  in 
Figure  11.  A  diagram  of  the  shock  tunnel  configuration  is  also  shown 
in  Figure  11.  It  is  clear  that  as  the  interface  diaphragm  location  was 
moved  closer  to  the  test  section,  the  helium  arrived  at  the  test  sec¬ 
tion  earlier. 

The  linear  relationship  between  the  ratio  of  helium  percentage, 
to  air  percentage  of  the  transient  gas  and  the  ratio  of  helium  intensity 
to  air  intensity  shown  in  Figure  10  indicates  that  the  interface  de¬ 
tection  system  is  indeed  measuring  the  per  cent  concentration  of  helium 
and  air.  At  a  helium  concentration  of  ten  per  cent  the  scatter  of  the 
data  in  Figure  10  indicates  an  accuracy  of  *4  per  cent  in  the  measurement 
of  helium  concentration.  The  data  in  Figure  11  show  the  ability  of 
the  system  to  detect  the  arrival  of  a  helium-air  interface. 

During  early  evaluation  of  the  system  the  radiation  detected  by 
the  helium  channel  during  a  tunnel  run  with  no  electron  beam  present 
was  investigated.  At  this  point  in  the  development  the  helium  channel 

O 

contained  a  filter  with  a  200  A  half  width.  It  was  found  that  on  some 
runs  the  helium  channel  exhibited  enough  background  radiation  to  cause 
questionable  results  for  low  percentages  of  helium.  The  background 
radiation  was  not  significantly  greater  than  the  emission  level  result¬ 
ing  from  low  percentages  (<  10  per  cent}  of  helium.  A  helium  filter 
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Fig.  11  Buffer  Gas-Driven  Gas  Interface  Arrival  Time  Study 
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with  50  A  half  width  was  placed  in  the  helium  channel  which  reduced  the 
background  radiation  to  an  acceptable  level.  The  data  in  Figure  11  were 

O 

taken  with  the  200  A  filter.  This  is  the  reason  for  the  difference  in 
intensity  ratio  for  air  in  the  early  part  of  the  trace.  This  also 
shows  the  need  for  careful  consideration  of  background  radiation  in  such 
systems. 

Another  result  of  system  imperfection  which  is  evident  in  both 
Figure  7,  page  36,  and  Figure  10,  page  41,  is  that  a  zero  helium  per 
cent  does  not  give  an  intensity  ratio  of  zero.  The  primary  reason  for 
this  is  that  the  helium  channel  indicates  some  output  due  to  the  much 
stronger  nitrogen  emission.  Also  the  imperfect  null  of  zero  frequency 
components  in  the  output  circuits  is  another  contributor. 
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CHAPTER  V 
CONCLUSIONS 


An  electron  beam  excitation  system  and  a  specialized  spectro¬ 
meter  for  measuring  transient  gas  concentrations  in  a  shock  tunnel  have 
been  developed.  Evaluation  data  demonstrated  that  the  system  is  capa-  ' 
ble-of  detecting  the  percentage  of  helium  in  a  transient  helium  air 
mixture  within  *4  per  cent  of  helium  percentage.  There  did  not  appear 
to  be  any  degrading  effects  due  to  variations  in  beam  current  (0-1  ma), 
pressure  (0.05-1  Torr),  or  temperature  (60-300°K).  Also,  the  time 
response  of  the  device  to  changes  in  mixture  concentration  indicates 
that  there  are  no  transient  effects  due  to  resonance  diffusion  in 
helium  or  collision  effects  between  molecules  of  the  same  gas  or  mole¬ 
cules  of  different  gases  in  a  helium-air  mixture  which  would  invalidate 
transient  concentration  measurements  taken  on  a  time  scale  of  a  few 
milliseconds. 

There  are  at  least  five  areas  which  warrant  further  study  if  the 
ultimate  potential  of  gas  concentration  measurements  with  the  electron 
beam  is  to  be  realized: 

1.  The  effects  on  electron  beam  emission  of  densities  equiva¬ 
lent  to  several  Torr  at  room  temperature  need  further  work.  There  are 
many  cases  where  density  and  gas  concentrations  could  be  of  value  if  a 
firm  upper  limit  for  density  were  established. 

2.  In  conjunction  with  high  density  studies  much  could  be 
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learned  about  the  resonance  diffusion  and  collision  effects  mentioned 
above  since  they  will  be  the  limiting  conditions  on  higher  density 
work. 

3.  It  is  evident  from  the  theoretical  calculations  of  electron 
beam  excited  emission  that  the  present  values  for  excitation  cross- 
section  and  even  lifetimes  will  have  to  be  improved  if  theoretical 
designs  of  such  systems  are  to  be  done  with  confidence. 

4.  A  need  which  is  too  often  put  aside  is  the  need  for  a  high 
current  (10  ma),  high  voltage  (50-100  KV)  electron  gun  which  is  easily 
adapted  to  experimental  apparatus.  High  quality  gun  operation  is  pre¬ 
requisite  to  any  viable  data  acquisition. 

5.  Additional  capabilities  in  the  instrumentation  system  would 
be  advantageous.  With  proper  interpretation  of  the  effects  of  beam 
spreading,  beam  attenuation,  and  catcher  cup  current  measurements, 
additional  dividers  could  be  added  to  take  the  ratio  of  absolute  inten 
sity  and  beam  current,  thereby  obtaining  the  density  of  helium  and  air 
directly.  Also,  a  beam  modulation  and  detection  system,  which  is 
discussed  in  the  Appendix,  would  reduce  the  effect  of  background 
radiation. 
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APPENDIX 

SIGNAL  TO  NOISE  IMPROVEMENT  IN  EXTERNALLY  EXCITED, 
MILLISECOND  DURATION,  RADIATION  PHENOMENA 


A  significant  improvement  in  accuracy  could  be  obtained  by 
eliminating  the  effects  of  background  radiation  from  the  electron  beam 
excited  radiation.  A  study  of  possible  methods  was  made  following  the 
work  discussed  as  the  subject  of  this  thesis. 

There  are  several  approaches  which  have  been  used  successfully 
to  improve  the  signal  to  noise  ratio  of  spectroscopic  data  taken  on  a 
one  second  or  greater  time  scale.  Quite  often  narrow  band  optical  and 
electrical  filters  are  sufficient.  However,  in  electron  beam  diagnos¬ 
tics  in  shock  tunnel  flows  a  significant  background  level  passes  through 
even  a  narrow  optical  filter  and  the  electrical  filter  must  have  a  band¬ 
width  sufficient  to  satisfy  millisecond  response  requirements. 

In  cases  where  the  noise  is  sufficiently  random  and  the  photon 
arrival  rate  is  low  enough,  single  photon  counting  techniques  have  been 
utilized  by  Arecchi,  Gatti,  and  Sona  (24).  However,  they  state  that 
for  counting  systems  with  microsecond  response  the  necessity  of  no 
appreciable  counting  losses  limits  single  photon  counting  techniques 
to  photoelectron  arrival  rates  less  than  10s  photoelectrons  per  second. 
Assuming  that  system  response  is  not  a  limiting  factor  it  is  revealing 
to  consider  the  number  of  photons  necessary  to  establish  a  valid  aver¬ 
age  data  point  every  100  microseconds.  Assuming  the  arrival  rate  of 
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photoelectrons  at  the  anode  has  a  Poisson  distribution  the  error  in 
determining  the  average  number  of  counts  for  a  fixed  sampling  period  is, 

o  »  ST.  (57) 

Where  o  is  the  standard  deviation  of  the  number  counted  per  sampling 
period  and  N  is  the  average  number  counted  per  sampling  period.  In 
order  to  obtain  a  standard  deviation  which  is  ten  per  cent  of  the  aver- 
age  number  counted  per  sampling  period,  a  total  of  10  photoelectrons 
must  be  counted  per  sampling  period.  For  a  typical  photocathode  quan- 

•X 

turn  efficiency  of  ten  per  cent  this  means  10  photons  must  be  incident 
on  the  photocathode  per  sample  period.  If  the  sample  period  is  taken 

7 

as  100  microseconds  then  the  required  photon  arrival  rate  is  10' 
photons  per  second.  The  photon  arrival  rates  inherent  in  the  data 
acquired  in  this  work  were  on  the  order  of  10^  photons  per  second. 
However,  while  the  information  rate  is  sufficient  the  unavailability 
of  submicrosecond  counting  equipment  precluded  immediate  investigation 
of  this  technique. 

Lock-in  or  synchronous  amplification,  where  the  excitation 
source  is  modulated  and  the  output  is  phase  detected,  has  been  used 
quite  successfully.  The  operation  of  synchronous  detection  is  related 
quite  succinctly  by  Schwartz  (25). 

"The  message  function  fm(t)  is  modulated  into  a  cosine  carrier 
at  the  transmitter  to  give 

fs(t)  =  fm(t)  cos  wQ  t.  (58) 


54 


AEDC-TR-6B-88 


At  the  receiver  the  incoming  signal  is  multiplied  by  a  cosine  refer¬ 
ence,  voltage  of  the  same  frequency  and  phase  as  the  carrier.  The 
result  is 


2  fmCt) 

fm(t)  cos'*  wQt  *  — 2~  (1  +  cos  2w0t).  (59) 

The  signal  is  then  passed  through  a  low -pass  filter  which  eliminates 
thedouble  carrier  frequency  term.  The  output  is  1/2  fm(t)." 

By  making  the  low-pass  filter  narrow  band  about  zero  frequency 
the  signal  to  noise  ratio  can  be  greatly  enhanced. 

■■  In  applications  of  interest  here  where  the  output  bandwidth 
cannot  be  made  arbitrarily  narrow  the  remaining  advantages  are  the 
3  db  increase  in  signal  to  noise  ratio,  as  indicated  by  Schwartz  (25), 
and  the  shifting  of  information  to  higher  frequencies  where  the  noise 
should  be  reduced  if  it  follows  the  common  1/f  energy  spectrum. 

Schwartz  (25)  also  indicates  that  the  greatest  advantage  to  synchronous 
detection  is  for  signal  to  noise  ratios  less  than  one. 

In  most  transient  cases  presently  of  interest  the  signal  to 
noise  ratio  is  greater  than  one.  Also,  a  brief  experiment  with  a 
Lock-In  amplifier  showed  that  the  phase  shift  through  the  electron 
beam  excitation  system  is  not  constant.  Therefore,  it  was  deemed 
prudent  to  investigate  a  less  complicated  amplitude  modulation  and 
detection  scheme  which  should  give  considerable  improvement  in  signal 
to  noise  ratio  since  it  incorporates  the  primary  advantage  of  synchro- 

r 

nous  detection  by  shifting  information  to  higher  frequencies. 
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A  block  diagram  of  the  signal  extraction  system  is  shown  in 
Figure  12,  For  convenience  a  45S  kilohertz  modulation  frequency  was 
selected  because  conventional  intermediate  frequency  components  are 
available.  The  modulation  signal  was  capacitance  coupled  into  the  con¬ 
trol  grid  of  the  electron  gun  through  a  455  kilohertz  intermediate 
frequency  transformer,  A  second  455  kilohertz  output  transformer  was 
used  at  the  output  to  tune  the  photomultiplier  output  circuit  to  the 
carrier  frequency.  Noise  was  introduced  by  modulating  a  neon  bulb 
biased  in  the  glow  regime  by  a  battery.  A  diode  detector  was  used  to 
provide  a  dc  signal  at  the  output. 

There  are  two  points  which  warrant  comment  regarding  this 
technique.  It  is  imperative  that  the  beam  be  modulated  cleanly.  Any 
amplitude  variation  of  the  carrier  at  the  input  will  be  transmitted 
through  to  the  output.  It  is  of  interest  to  note  that  when  a  diode  is 
used  as  a  detector  and  the  carrier  is  greater  in  magnitude  than  the 
noise  frequencies  near  the  carrier  frequency,  the  carrier  switches  the 
diode  and  therefore  tends  to  enhance  only  information  in  phase  with  the 
carrier.  For  this  condition  amplitude  detection  is  somewhat  analogous 
to  synchronous  detection. 

Using  the  system  shown  in  Figure  i2,  signal  co  noise  improvement 

3 

of  10  was  readily  obtained.  A  plot  of  signal  ro  noise  improvement 
versus  noise  frequency  is  shown  in  Figure  13,  page  57, 


56 


AEDC.tr.  68-88 


Fig.  12  Amplitude  Modulation  and  Detection  System 
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